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Abstract 

We investigate what happens to the third order ferromagnetic phase transition displayed by the 
Ising model on various dynamical planar lattices (ie coupled to 2D quantum gravity) when we intro- 
duce annealed bond disorder in the form of either antiferromagnetic couplings or null couplings. We 
also look at the effect of such disordering for the Ising model on general (f>^ and (f>'^ Feynman diagrams. 
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In an earlier paper we observed that the Ising antiferromagnet would not undergo a phase tran- 
sition on dynamical 4>^ or 0"* planar graphs or on dynamical planar triangulations, essentially because 
of frustration, but would display a transition on dynamical planar quadrangulations. Here we address 
a related question, namely: what degree of annealed bond disorder is necessary in order to destroy the 
ferromagnetic phase transition on a dynamical lattice? For simplicity we consider models in which there 
are two possible bond factors, + J appearing with probability p and either — J or appearing with prob- 
ability {1 — p). We consider primarily planar (jr" and (f/^ graphs, along with planar triangulations and 
quadrangulations but also look more briefly at general (fi^ and 0'* graphs of arbitrary topology. We are 
thus further complicating the annealed connectivity disorder of the Ising model coupled to 2D gravity 
(ie living on a planar dynamical graph or triangulation) by introducing an additional annealed bond 
disorder. Such annealed bond disorder has been considered some years ago for the case of fixed regular 
lattices where it was shown that the critical exponents underwent a Fisher renormalization |3| for 
models with an initial ^ a > 0. For the two dimensional Ising model it was shown that introducing a 
finite fraction of either antiferromagnetic or null bonds suppressed the initial phase transition. For the 
Ising model on the dynamical lattices we consider here with a = —1 the considerations in will not 
apply, and it is not immediately obvious what behaviour to expect. 

To investigate the problem we consider the following Ising partition function on an ensemble of planar 
random graphs G" with n vertices. 

ZniPJ) = E E / ^(•^)'^-^ f E '^^'^j I (1) 

G" {a} \ <ij> ) 

where is the connectivity matrix for a given graph and P{J) is the annealed bond distribution. For 
the case of a trivial P{J) = 5{J — 1) summing over the number of vertices gives 

^(^'5) = E((r3^j Zr^m (2) 

where c — exp(— 2/3), so equ.(|l|) reduces to the Ising model solved by Boulatov and Kazakov in Q by 
noting the equivalence of equ.(ph to the free energy of a two matrix model 



(3) 



where a, 6 = ± and (f>^ are N x N matrices. Borrowing the notation of ||] the inverse propagator K^^^ 

is given by 

where the coupling g = exp(2/3) = 1/c and the potential is 

y(0+,0_) = ^(0^-i-03_) (5) 

for (j)'^ graphs and 
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= + (6) 

for (j)* graphs. In the 0^ case, after some rescalings, the matrix model action can be written as 



^i.e. In the original ordered version of the model 
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and the matrix integrals in the partition function can be carried out using the methods of Q to obtain 
an exact expression for the partition function, even in the presence of an external field. The (non- 
Onsager) critical exponents that are derived from this expression are consistent with those obtained 
in the continuum formulation using either lightcone or conformal gauge Liouville theory Q]. In the 
second paper of Q the model was also solved on a cj)^ graphs and the same exponents obtained, thus 
confirming universality. This dynamical lattice universality is further supported by the results of |^ for 

graphs without tadpoles and self-energy diagrams, which again give the same set of critical exponents. 
A considerable body of numerical evidence supporting these results has been accumulated to date. 

Let us now consider the effect of introducing a non-trivial bond distribution P{J) in the above 
calculations. If wc want to examine the effect of introducing ferromagnetic bonds with a probability p 
and antiferromagnctic bonds with a probability (1— p) we should take P{J) = p 6{J —l) + {l — p) 5{J +1). 
This has the effect of modifying the inverse propagator to 

K-j- w r ^ ) t f ) (8) 




V5 ^ / 

We can now invert the propagator |^ and rescale the resulting action to a form resembling that of equ.(|^) 

S^cl?, + 4>l~2c{p)4>+4>_-^{^X + ct^t) (9) 



where 



and 



C{P) = ^+^7"^^ (10) 



P .9 + (1 - P) 



X = 2y^—±l^fi-^, (11) 



4( ff^-l)2(2p-l )^A 

g{p9 + 1 -py 

As we have not changed the form of the action with these manipulations, merely the interpretation of c 
and A which are parameters in the action, the equations for determining the critical point arc unchanged. 
We thus have, substituting c{p) for c in the appropriate equations in [Q, 

P+ (1 -P) 9c ^ 1 .-^2) 
pgc + {l-p) 4 ^ ' 

where gc is the critical value of exp 2(3. Solving this wc find 

exp2/3, = ^^ (13) 
bp — A 

which tells us that (3c ^ oo as p ^ 0.8. In other words the critical temperature of the theory goes to 
zero as the probability of ferromagnetic bonds drops to 0.8 at which point the phase transition vanishes, 
leaving only a disordered phase. The value of p at which this happens is non-universal as can be seen by 
considering 0^ lattices. Again the form of the action is unchanged by introducing the bond disorder and 
c in the original ferromagnetic theory is replaced by c(p) to give the equations 

P+a-p) gc ^ a/28-1 

pgc + {l^p) 27 ^ ' 

for a pure <p'^ lattice, which gives the vanishing point p ~ 0.86 and 

p + {1 - p) gc _ 23 



pgc + {l-p) 108 



(15) 



for a (/)^ lattice with no tadpoles and self energies, which gives the vanishing point p ~ 0.82. It is perhaps 
worth remarking that these numbers are larger than the percolation thresholds on the appropriate graphs 

^This is singular for p = 1/2, but we run into trouble before reaching this probability. 
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without matter where they have been calculated (p ~ 0.78 on the (jy^ graphs with tadpoles and self-energies 
and 2/3 on the graphs jl^), though the back reaction of the matter might be expected to modify the 
geometry and hence these values at the Ising critical points. 

We can play similar games in looking at the effects of bond disorder in bonds of the same sign in the 
model. If we take, for instance, P{J) — p 5{J — 1) + (1 — p) 5{J — 2) the propagator becomes 

K-a,'- P( f +i^-P)( I I) (16) 



^7^ V9 J \-a 9 

an we can absorb the change in a rescaling of A and a c{p) given by 

On the lattice for instance the critical value of g will thus interpolate smoothly between 4 and 2 as 
p changes from 1 to 0. For P{J) of the form p 5{J — 1) + (1 — p) 6{J — A) the interpolation would be 
between 4 and 4^^/^^. 

The case of bond dilution, in which null bonds are inserted with probability {1 — p) is particularly 
intriguing. We have P{J) — p S{J — 1) + (1 — p) S{J) in this case, so 



f( ^ ^1 +('-")( 1 1 I (18) 



which in turn leads to a c{p) of the form 

r(r>) — 

pV9+ (1 -P) 



. ^ p/yg + i^-p) .,„^ 
c{p) = — _ , ^ (19) 



and 

3^ ^ 4p^(ff-l)^(p(.g-H) + 2(l-p)Vg)^A ^^q) 
If we now solve the equation for the critical point c(p) = 1/4 on (j)^ graphs, we find that 



exp(/3.) . 3(l-p)W9(l-p)-+i6p^ (21) 
2p 

which gives a /3c — *■ oo only as p — > 0. A similar result applies for the bond dilute model on (jr' graphs. 
The critical temperature goes to zero and the transition vanishes only when the probability of non-null 
bonds goes to zero. This is in direct contrast to the case of annealed disorder with dilute bonds on a fixed 
lattice where the transition vanishes for a finite value of p that is lower than in the case of competing 
bond signs and close to the percolation threshold for the lattices It is rather surprising that the 
transition should persist on a dynamical lattice down to zero active bond concentration - perhaps the 
dynamical nature of the lattice allows sufficient communication between the active bonds to generate a 
transition by some sort of clustering mechanism. In the fixed lattice case small, but definite, correlations 
are present between like bonds, so such an effect could well be magnified on any sort of dynamical lattice. 

In 1^ we observed that it was also possible to write down directly matrix models for Ising spins on 
dynamical triangulations 

S^-tr(- — rrir. T-^' + T^ — rrir. rD^ + sys + SoA (22) 

g V2cosh(/3)(l c*) 2 cosh(/3)(l - c*) ' / ^ ' 

where S represented edges of triangles with the same spins at each end and D edges of triangles with 
different spins. The coupling c* was the "dual" to c, (1 — c)/(l -|- c). The model could be rescaled 
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into the 0(1) representation of the Ising model written down in [pl]| and gave the correct dual critical 
temperature. We can generalize the arguments of to an action of the form 

S ^ —tr ( -A{(3)S^ + -B{p)D^ + /?, + SdA (23) 
5 V2 2 J 

which, using the results of lITTI , will display Ising critical behaviour if A and B satisfy 

7B{l3f + &A{l3f ~ UA{(i)B{l3) = 0. (24) 

The propagator is diagonal in such models, which makes it particularly easy to write down the effect of 
non-trivial P(J)s. For P( J) ^ p S{J -l) + {l~p) S{J+1) we have lM(/3,_p) = cosh(/3)(l + (2p-l)c*) and 
1/B{(3,p) = cosh(/3)(l— (2p— l)c*). Dividing through by A^ in equ. (p4|) we see that Ising critical behaviour 
is possible if the ratio A/B takes the same value as in the case p = 1, namely (14 - V7)/(13 + \/7). We 
thus find that 

e* = ^^^ (25) 
2p - 1 27 ^ ^ 

which gives a minimum value of p, since we know that a ferromagnetic transition must have c* < 1. The 
transition is pushed to zero temperature (c* = 1) at p ~ 0.58, leaving only a disordered phase. This is 
entirely analogous with the model on cf)^ and 0^ graphs. Taking the ratio of the couplings in the bond 
dilute case where 1/A{l3,p) = pcosh(/3)(l + c*) + {I - p) and l/P(/3,p) = pcosh(/3)(l - c*) + (1 - p), 
gives the following equation for the critical temperature 

cosh(/3c) - 6.291 sinh(/3c) ~ (p - l)/p. (26) 

This has solutions for arbitrarily small p so, just as for the dilute bond Ising model on planar cl)^ and 
graphs, the phase transition only vanishes as p ^ 0. 

The matrix model for dynamical quadrangulations considered in , which can be written in a rescaled 
form as 

S = r775Tn T^' + Z r775Tn T^' + 1^^ + 7^' + U^DSD + 2S^D^)) , (27) 

g V2cosh(/3)(l +c*) 2 cosh(/3)(l - c*) 4 4 2 jj^ \ j 

can also be modified in a similar fashion to the triangulation model to include bond disorder: 

S^—tr( -A{I3)S^ + -B{P)D^ + -S^ + -D^ + -{SDSD + 2S^D^)) . (28) 
g \2 2 4 4 2 / 

If we include antiferromagnetic bonds, P{J) = p 5{J + 1) + (1 — p) &{J — 1), and demand that the ratio 
B/Ahe the same as in the original model we find 

c* = (29) 
2p - 1 4 ^ ' 

Noting that c* should be less than 1 again gives a critical value of p = 5/8 for the vanishing of the 
ferromagnetic transition. The dynamical quadrangulation model will also have an antiferromagnetic 
transition, which will appear for p < 3/8 by similar arguments upon replacing c* by — c* (ie j3 by —(3). 
For 3/8 < p < 5/8 there would appear to be no transition at all. The bond dilute model gives the 
following equation for the critical temperature 

cosh(/3,) - 4sinh(/3,) = (p - l)/p (30) 

so the critical temperature will only go to zero as p — > 0. 

The planar (j)^ and 0'* graphs and their duals that we have looked at here can be considered as an 
ensemble of "fat" or ribbon graphs generated by the perturbative N x N matrix integrals in the limit 
A'^ — > oo, with the spin models living on them. In a recent paper |^ the opposite limit of iV — > 1 was 
considered, which gives an ensemble of standard "thin" Fcynman diagrams. In an annealed ensemble of 



4 



such graphs, which are locally tree like, the Ising model displays a mean field transition at the Bethe 
lattice values of Qc, namely 3 for cj)'^ graphs and 2 for graphs. Introducing a coupling distribution of 
the form P{J) = p 6{J — 1) + (1 — p) S{J + 1) on thin (j>^ graphs and scaling the action the form 

S = l{cf>l + cl^l) - c{p)cl,+c^- '^-) (31) 

gives the following saddle point values for , 0_ 



(1 -P) 



, , 1 + g ± v/(g + l)(g(4p - 3) - (4p^ 

^ 2((l-p)+pg) (^2) 

(along with a trivial zero solution and a solution with (f>^ and 0_ interchanged) ||] . The low temperature 
ordered phase is given by the second solution with differing , 0- , so this is only possible when the 
expression inside the square root is positive. The critical gc is that for which the square root is zero, 
namely 

9c = (33) 

Once again the critical temperature goes to zero at a finite value of p, leaving only the disordered phase 
for p < 3/4. The rather peculiar behaviour of the dilute model with P(J) = p S{J — 1) + {1 — p) S{J) 
also appears to persist on thin graphs, where we have the following saddle point solutions for 0+, ^_ 

^ . (3-l)P 



(1 -p)\/5+W 



2Vg(l -P)+ Pjg + 1) ± Vp'g'' - 2(3p^ - 4p + 2)g - 8p(l - p) yg - 3p^ 

2(Vg(i-p) + .gp) ^ ^ 



In this case we find that it is possible to obtain a positive term inside the square root for all values of p, 
with the critical temperature going to zero only as p ^ 0. 

The saddle point equations for general 0"* graphs may be derived in a similar fashion from the action 

S = \{cf>l + cl^l) - c(p)</)+0- - + '^-) (35) 

for the various c{p). For purely ferromagnetic interactions, c(p) = 1/g, the high and low temperature 
saddle point values are 



271 ' 



1 - yi ^ 4/g2 

71 



(36) 



We find that for antiferromagnetic bonds, P{J) = p S{J — 1) + (1 — p) S{J + 1), the ordered phase is only 
possible for p > 2/3, with the critical value of g being given by 



3p- 1 
3p- 2' 



(37) 



Bond dilution gives a transition that persists for all finite p in this case too. 

In summary, we have investigated the effects of including annealed bond disorder and dilution on the 
Ising model in ensembles of planar cj)^ and 0^ graphs, triangulations, quadrangulations and general (j)'^ and 
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(j)'^ graphs. We find that a fixed finite probability of antiferromagnetic bonds is sufficient to suppress the 
Ising transition in aU cases, but that the transition persists down to zero active bond probabihty in the 
case of bond dilution, with the transition temperature going to zero as p — > 0. It would be interesting to 
attempt to understand the persistence of the bond dilute transition by looking at bond/bond correlations 
in the various models in a similar fashion to the fixed lattice work in It might also be possible to 
extend the calculations to bond distributions other than the sums of delta functions considered here. 
From the numerical point of view it would be an amusing exercise to perform simulations to attempt 
to verify the predictions for the critical concentrations with antiferromagnetic bonds and to see if the 
transitions in the bond dilute models really are as tenacious as the calculations here suggest. 

We have not considered quenched connectivity or bond disorder at all in this paper. There is no 
problem in dealing with this on the thin graphs and the calculations in |^ indicate that interesting 
effects such as a spin glass phase may appear on the introduction of quenched bond disorder - or simply 
antiferromagnetic couplings. For "fat" graphs we immediately run into problems with the c = 1 barrier 
when we employ the replica trick to calculate logarithms of partition functions. For the Ising model we 
can only caculate with n = 1, 2 replicas, for instance ^ In view of the importance of the replica trick in 
this domain any insights on matrix models for c > 1 are likely to have some bearing on the treatment of 
quenched disorder. 

D. Johnston was supported at Orsay by an EEC Human Capital and Mobility Fellowship and an 
Alliance grant. 
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